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Abstract 

This article describes the fluid dynamics video, "A Numerical Sim- 
i ulation of a Plunging Breaking Wave" , which was submitted to the 

gallery of fluid motion at the 2009 APS/DFD conference. The simula- 
tion was of a deep-water plunging breaking wave. It was a two-phase 
calculation which used a Volume of Fluid (VOF) method to simulate 
the interface between the two immiscible fluids. Surface tension and 
If} viscous effects were not considered. The initial wave was generated by 

applying a spatio-temporal pressure forcing on the free surface. The 
video shows the 50% isocontour of the volume fraction from several 
t-H different perspectives. Significant air entrainmcnt is observed as well 

as the presence of stream- wise vortex structures. 
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The computational domain, shown schematically in Fig. [I] moves with the 
linear crest speed of the wave, U, and has dimensions [2ir, tt/2, 2ir], with 
[1024,256,512] grid points. Periodic boundary conditions are used in the 
horizontal directions and free-slip boundary conditions in the vertical. The 
Froude number, Fr = U(kfg) 1 / 2 , is unity, where g is the acceleration due to 
gravity and k is the wave number. The density ratio between the two fluids is 
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1 Formulation 
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Figure 1: Schematic of computational domain, g = [0,0, —g]. 



1000:1. The total distance traveled by the wave, X in the laboratory frame, 
is equivalent to Ut in the simulation. Gravity acts in the —z direction. 



2 Numerical Method 

Let u* denote the three-dimensional velocity field, p the density, and p the 
pressure, all functions of space (x*) and time (f), where an asterisk denotes 
a dimensional quantity. Time, space, velocity, density and pressure are 
normalized as follows: 

t = t*Uk , Xi = x*k , Ui = y- , p=— , v=^-j^k. (1) 

U p w p w U z 

where p w is the density of water, k is the wave number, and U = (g/k) 1 ^ 2 . 

In the Volume of Fluid method (Rider et al. (1994)) the fraction of fluid 
that is inside a cell is denoted by <fi. By definition, = for a cell that is 
totally filled with air, and 0=1 for a cell that is totally filled with water. 
The density expressed in terms of 4> is 

p = (l_A)0 + (2) 

where A is the density ratio between air and water. 
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The non-dimensional governing equations are: 
Momentum: 



duj 

at 



d {u k Uj 
dx k 



VOF: 



1 dp 

pdxi 



dt Uk dxk 



0. 



(3) 



(4) 



The relevant non-dimensional parameter is the surface Froude number, Fr = 
U (k/g) 1 / 2 . The temporal integration is handled with an explicit RK2 scheme, 
and the advective terms with the flux based limited QUICK scheme of 
Leonard (1997). The VOF algorithm uses the operator-split method of 
Puckett et al. (1997). As discussed in Dommermuth et al. (1998) the di- 
vergence of the momentum equation, Eq. Q, combined with the solenoidal 
constraint dui/dxi = 0, provide a Poisson equation for the dynamic pres- 
sure. 

The last term on the r.h.s. of Eq. (|3]) is the atmospheric pressure forcing 
term. In the simulation discussed here the following forcing function was 
used for t < 4ir 
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[Aqcos (a?) + A r 



1 



cos 



where 
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dVA r (an) 



(5) 



(6) 



Here, Aq = 0.02 and A r is a uniform random disturbance that has been 
passed through a low-pass filter. 



3 Flow description 

Fig. [2] is a plot the the total energy, 

E(t) = f dVp(t) {U t {t)U t {t) + gz) , (7) 

over time. Four distinct stages are evident in Fig. [2j they are: A. Atmo- 
spheric forcing; B. Potential flow before breaking; C. Breaking which con- 
sists of plunging, spilling and splash-up events; and D. Potential flow after 
breaking. 
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Ill 
A: Atmospheric forcing B: Potential flow before breaking 
C: Plunging & Spilling D: Potential flow after breaking 
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Figure 2: Total Energy, J v dVp (UiUi + gz) 

4 Video 

The full size video, mpeg2 encoded, is approximately 75Mb ( [download] ) . 
The web size video, mpegl encoded is approximately WMb (download). 



The videos are also available at www . saic . com/ nfa} 



Scenes 1 and 2 in the video start near the end of stage A (see Fig. [2J and 
show the evolution through stages B and C and end early in stage D. Scene 3 
starts after the initial plunging event, early in stage C and ends late in stage 
C. Scene 4 starts near the end of stage B, and ends in the middle of stage 
C. Scene 5 (Vortex Tubes) starts near the end of stage B and pauses/ends 
early in stage C. In scenes 1, 2 and 3 the 50% isocontour of the volume is 
shown. In scenes 4 and 5 the magnitude of velocity projected onto the free 
surface is shown in selected regions. 

SCENES: 

1. Tank View- A perspective view with walls added to provide context 
to the breaking events. After the initial plunging event, a jet is formed 
(splash-up). After the secondary jet impinges on the free surface the 
remaining portion of stage C is a sequence of weaker and weaker plung- 
ing and spilling events, until the flow returns to a potential flow. 



2. 4-way View- Here the evolution is viewed from four angles simultane- 
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Figure 3: Vortex tubes on the underside of the initial oval of air entrained 
by the first plunging event. 

ously. Note the large pocket of air that is entrained (lower right view) , 
and the irregularity on the front face of the wave (lower left view) . 



3. Side (bottom) Zoom- This scene is a magnified view of the lower right 
frame from Scene 2, and shows the formation of the jet during splash- 
up. Significant air entrainment is observed shortly after the impact 
of the splash-up jet. This air entrained by the splash- up jet persists 
for a longer time than the initial ovular pocket of air entrained by the 
initial plunging event. 



4. Velocity Magnitude- This is Scene 4 with the magnitude of the velocity 
projected onto the free surface. 

5. Vortex Tubes- Scene 4 paused to show the vortex tubes, and hairpin 
like structures on the back of the initial oval of air entrained by the 
plunging event. Shown here as Fig. [3] 
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